Genome-wide detection of human variants that disrupt intronic branchpoints.

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Trans-ethnic genome-wide association study of severe COVID-19.

COVID-19 has caused numerous infections with diverse clinical symptoms. To identify human genetic variants contributing to the clinical development of COVID-19, we genotyped 1457 (598/859 with severe/mild symptoms) and sequenced 1141 (severe/mild: 474/667) patients of Chinese ancestry. We further incorporated 1401 genotyped and 948 sequenced ancestry-matched population controls, and tested genome-wide association on 1072 severe cases versus 3875 mild or population controls, followed by trans-ethnic meta-analysis with summary statistics of 3199 hospitalized cases and 897,488 population controls from the COVID-19 Host Genetics Initiative. We identified three significant signals outside the well-established 3p21.31 locus: an intronic variant in FOXP4-AS1 (rs1853837, odds ratio OR = 1.28, P = 2.51 × 10-10, allele frequencies in Chinese/European AF = 0.345/0.105), a frameshift insertion in ABO (rs8176719, OR = 1.19, P = 8.98 × 10-9, AF = 0.422/0.395) and a Chinese-specific intronic variant in MEF2B (rs74490654, OR = 8.73, P = 1.22 × 10-8, AF = 0.004/0). These findings highlight an important role of the adaptive immunity and the ABO blood-group system in protection from developing severe COVID-19.

SARS-CoV-2 orf1b Gene Sequence in the NTNG1 Gene on Human Chromosome 1.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-sense single-stranded RNA virus. It is contagious in humans and is the cause of the coronavirus disease 2019 (COVID-19) pandemic. In the current analysis, we searched for SARS-CoV-2 sequences within the human genome. To compare the SARS-CoV-2 genome to the human genome, we used the blast-like alignment tool (BLAT) of the University of California, Santa Cruz Genome Browser. BLAT can align a user sequence of 25 bases or more to the genome. BLAT search results revealed a 117-base pair SARS-CoV-2 sequence in the human genome with 94.6% identity. The sequence was in chromosome 1p within an intronic region of the netrin G1 (NTNG1) gene. The sequence matched a sequence in the SARS-CoV-2 orf1b (open reading frames) gene. The SARS-CoV-2 human sequence lies within non-structural proteins 14 and 15 (NSP14 and NSP15), and is quite close to the viral spike sequence, separated only by NSP16, a 904-base pair sequence. The mechanism for SARS-CoV-2 infection is the binding of the virus spike protein to the membrane-bound form of angiotensin-converting enzyme 2 and internalization of the complex by the host cell. It is probably no accident that a sequence from the SARS-CoV-2 orf1b gene is found in the human NTNG1 gene, implicated in schizophrenia, and that haloperidol, used to treat schizophrenia, may also be a treatment for COVID-19. We suggest, therefore, that it is important to investigate other haloperidol analogs. Among them are benperidol, bromperidol, bromperidol decanoate, droperidol, seperidol hydrochloride, and trifluperidol. These analogs might be valuable in the treatment of COVID-19 and other coronavirus infections.